"0","all.sims <- rbind(lowE,intermediateE,highE)

limits <- aes(ymax = EsN + Es.se,ymin=EsN - Es.se)
p <- ggplot(all.sims,aes(N,EsN,colour=model))
p  + geom_line() + geom_point() + geom_errorbar(limits,width=0)"
